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ABSTRACT 


An exploratory study has been made of the stage separation of 
parallel staged reusable launch vehicles. Static longitudinal aerody- 
namic data were obtained for both stages of a representative tvo-stage 
rocket-powered reusable launch vehicle concept when the stages were in 
close proximity to each other. The effects of vertical spacing, 
longitudinal spacing, and incidence angle were determined at Mach 
numbers of 3 and 6. Using the wind-tunnel data together with estimates 
of the dynamic derivatives obtained for similar type vehicles and also 
the vehicle and initial trajectory input data, the implication of the 
interferences on the individual vehicle dynamics was obtained at Mach 
numbers of 3 and 6 by numerically integrating the equations of motion 
for both vehicle stages. 



CHAPTER I 


INTRODUCTION 

During the past several years, interest in the potential of the 
recovery and reuse of launch vehicles has resulted in a number of studies 
of various types of recoverable systems. The spectrum of concepts has 
varied from recoverable ballistic to winged reusable air-breathing 
vehicles . The latter type and most of those which employ aerodynamic 
lift during the ascent trajectory would probably utilize multiple stages 
which are arranged in parallel for launch. If staging must occur within 
the sensible atmosphere, aerodynamic interferences between the two 
vehicles may result in significant effects on the behavior of each 
vehicle during the staging maneuver. The author has therefore undertaken 
an exploratory investigation to ascertain these interferences and to 
interpret their influence on the staging maneuver. The present paper is 
an initial attempt at exploring the complex problem of separating 
parallel arranged stages of a reusable launch vehicle system. 

In order to provide meaningful information on the magnitude and 
character of the staging problem it is necessary to obtain aerodynamic 
test data on each of the two vehicles at conditions which might occur 
during staging and then to estimate the effects of these forces and 
moments on the relative motion of the two vehicles. For the present 
investigation, which is primarily exploratory, static longitudinal data 
were obtained on two representative stages of a reusable launch vehicle 
system for a range of vehicle positions and attitudes to obtain some 
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understanding of the aerodynamic interference phenomena. These data 
were then used in numerically integrating the equations of motion of 
both vehicles to determine the relative behavior of both stages during 
a staging maneuver. 

The experimental aerodynamic data were obtained for a model of a 
representative rocket -powered two-stage launch vehicle system. The first 
stage was a simplified wind-body configuration combination, whereas the 
second stage was a lifting body. The measured data consisted of static 
longitudinal aerodynamic forces and moments which were obtained simul- 
taneously for each stage with strain-gage balances mounted on separate 
stings. Vertical and longitudinal spacing, as well as incidence angle, 
were varied for a range of angles of attack at Mach numbers of 3 and 6. 
The test results are used to illustrate the effect of these variables 
on the longitudinal stability of the two stages when in close proximity 
to each other. 

The analysis of the staging maneuver was made using the coupled 
linearized longitudinal equations of motion for both stages in order to 
calculate their relative behavior and the influence of aerodynamic 
interference. In these equations the general load terms were replaced 
by the experimental aerodynamic data and by estimates of the damping 
derivatives obtained fcr similar type vehicles at interference-free 
conditions. For selected values of stage characteristics and initial 
conditions the equations were solved numerically. Typical results are 
presented to show seme of the effects of vehicle parameters as they 
influence the safe separation or collision of the two stages. 



CHAPTER II 


APPARATUS AND TESTS 
Tunnel 

The wind-tunnel investigation was conducted in the 2-foot hyper- 
sonic facility at the National Aeronautics and Space Administration's 
Langley Research Center. The 2-foot hypersonic facility, described in 
reference 1, is a continuous, closed-circuit ejector type wind tunnel. 
The tunnel has a 24-inch by 24-inch by inch test section and provides 
for testing in the Mach number range of 3 to 6 at approximate Reynolds 
numbers between 2.0 x 10^ and 0.5 X 10^. The air supply for the tunnel 
is obtained from two 100,000 CEM compressors connected in parallel and 
operated at a pressure ratio of 4, and a 10,000 CFM auxiliary compressor 
feeding the main nozzle of the tunnel. The stagnation temperature of 
the tunnel is controlled by a 1,000-horsepower nichrome -tube -type heater 
located just upstream of the nozzle. A schematic drawing of the main 
test section is presented in Figure 1. The two-dimensional nozzle, 
diffuser, and ejector are individually supported on screw jacks and can 
be moved in the vertical direction. Thus, adjustments are available to 
set the desired tunnel conditions by proper positioning of the various 
tunnel components. A silica gel drying system is located in the circuit 
to monitor the flow into the test section at a dewpoint of -15° F or 
less, which is adequate to insure no moisture condensation shocks which 
would adversely affect the consistency of the te6t results . 
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Models 

Details of the launch vehicle configuration selected for the 
present study are shown in Figure 2. The launch vehicle consisted of a 
simplified wing-body first stage with a lifting-body second stage. The 
second stage was placed on top of the first stage, with the flat bottom 
parallel with the first-stage wing upper surface. The longitudinal 
location of the second stage was such that its moment reference center 
was behind the first-stage moment reference center in the normal carrying 

A v 

position, that is, = -0.051. Photographs of the two-stage system 

f l 

prior to staging and at staging conditions are shown in Figure 5 . 

The first stage consisted of a semicylindrical fuselage with an 
ogival forebody and delta wing. The wing had 65° leading-edge sweep 
while the airfoil section was a half -diamond section that had a maximum 
thickness of U percent of the local chord at the ^t-0 -percent chord 
station. The wing was flat on the upper surface, thereby giving the 
wing negative camber. No longitudinal or vertical control surfaces were 
provided on first stage for the present tests. 

The second-stage vehicle was a lifting body with cross-sectional 
shape progressions along the X-axis shown in Figure 2c . The second stage 
had essentially a flat bottom and was provided with vertical stabilizing 
surfaces . Pitch controls were provided on the second stage which 
consisted of upper and lower surface flaps located near the base of the 
body (see Fig. 2c). The ratio of the planform area of the pitch controls 
to the total planform area was 0.115 for the upper surface controls and 
0.187 for the lower surface controls. 
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Support Mechanism 

Separate sting supports were provided for the first and second 
stage, with the vertical movement between the stages being provided oy 
the support system to which the stings were attached (see Fig. 4). 
Longitudinal movement and incidence angle between the stages were 
provided using spacers and sting adapters on the upper sting support. 

The complete support apparatus was attached to an arc strut which varied 
the angle of attack of the vehicles. 

Tests 

The wind-tunnel tests were conducted at nominal Mach numbers of 
3.0 and 6.0. Position variables Ax and Az, as well as the relative 
incidence angle, i, were varied for an angie-of -attack range of approxi- 
mately -9° to 12° (see Fig. 5). The geometric variables Ax, Az, and 
i were systematically varied to obtain the aerodynamic data needed as 
input to the two-body trajectory simulation computer program. 

Static aerodynamic force and moment data were simultaneously 
obtained for the first and second stages by use of individual internal 
six-component strain-gage balances. No composite configurations, that 
is, with the first stage and second stage connected, were tested. 

All data were obtained with the model smooth; that is, no boundary- 
layer transition strips were used. At the Reynolds numbers of these 
tests laminar flow may be expected to exist over almost the entire length 
of the models. Individual vehicle angles of attack were corrected for 
balance and sting deflection under load. No base drag corrections were 
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made for either stage. The average test conditions and Reynolds number 
variations were as follows : 


Mach 

Stagnation 

Stagnation 

Reynolds number 

Number 

pressure 

temperature 

per foot 

3 

0.5 atm 

560° R 

0.8 x 10 6 

6 

3.0 atm 

76o° R 

1.0 x 10 6 



CHAPTER III 


METHOD OF ANALYSIS 

The separation maneuver of two parallel stages is illustrated 
schematically in Figure 6. At release, the second stage would be at an 
initial spacing distance and attitude with respect to the first stage. 

A trapeze or similar mechanism could be employed to achieve the desired 
release conditions . Potential problems that may arise during the staging 
sequence are illustrated by the lower two sketches of Figure 6. The 
divergence of the center of gravities alone would not imply safe separa- 
tion since the second stage may rotate into the first stage. Conse- 
quently, realistic analysis can only be accomplished if the separation 
maneuver is regarded as the motion of two finite rigid bodies . 

Two parallel lifting stages, separating from each other aerody- 
namically, represent a complex dynamic system. For the complete 
description of this system, not less than 12 degrees of freedom must be 
taken into account. While this is possible, at least theoretically, 
the final accuracy of the analysis will depend on the precise analytical 
representation of the aerodynamic forces and moments. Accordingly, the 
general motion of the dynamic system must be restricted to cases for 
which the approximate behavior of the static aerodynamic coefficients 
is presently known or can easily be determined experimentally. 

It was assumed that the controls of both vehicle stages were 
fixed and that the separation maneuver takes place above a flat non- 
rotating earth . Consequently, the equations needed to study the 
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motion of either stage are those which are derived in the Appendix, 
equations (A-27), (A-28), (A-29), and (A-30) . These equations can be 
simplified considerably in the case of longitudinal (symmetric) motion 
where it is assumed that the launch vehicle is flying with no sideslip 
and with wings level and that separation of the two stages takes place 
in the common plane of symmetry. With these restrictions, the equations 
necessary to study the motion of either stage become 

m(u + wq) = F* - mg sin 0 
m(w - uq) - F' + mg cos 0 

I yy < ^ 

0 = q 

x' = u cos 0 + w sin 0 
z ' = -u sin 0 + w cos 0 

The above set of six equations is completely general for symmetric 
vehicles and valid even for large disturbances. However, in the develop- 
ment of the expressions for the force and moment terms, the following 
assumptions are made: 

1. During separation the effects of speed changes on the aero- 
dynamic force coefficients are neglected. 

2. The axial -force coefficient of either stage is not affected 
by the proximity of the other stage. 

3. Cross coupling between the two stages takes place only through 
the normal -force and pitching -moment coefficients. 
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With the above assumptions, the force and moment terms may be 
expressed as 


Pi = | pv/sCGj, COS € - C A ) 

F> = | pV c 2 s(c T sin e - C N + ^ C Nq q) (2) 

My = | pV c 2 Sc(^ C T + C* + gSr C^q + C^w) 


In the above equations the rotational stability derivatives 
Cjj and C^, and the derivative C^. were assumed to be constants 
for each stage during the separation maneuver. However, the influence 
of these derivatives on the separation maneuver can be investigated by 
varying the magnitude of these quantities in the numerical analysis. 
Orientation of the thrust vector is specified by e, the angle between 
the X-axis and the thrust vector, and d^,, the perpendicular distance 
from the center- of gravity to the thrust line. These quantities, as 
well as the thrust coefficient Qj,, are constant for each stage. The 
static aerodynamic coefficients, C A , Cjj, and represent the 

experimental data measured on each stage and are defined separately 


as follows: 
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C A,1 " C A,l( a l) 

C A,2 = C A,2^ a 2^ 

= Cjj^(<x-^,<X2,Ax,Az) 
C N,2 = 

c m,l = Cm^Cai^AXjAz) 
C m,2 = G m,2(a 1 >oL 2 ,ZSx jr ^z) 


The separation variables a^, ag, Ax, and Az are illustrated in 
Figure 5. 

The equations of motion were integrated numerically on a digital 
computer. The values of the aerodynamic coefficients were obtained by 
linear interpolation between the discrete points of the tabular functions 
at each integration step. The atmospheric density, p, in equation (2) 
corresponds to that given by the U.S. Standard Atmosphere 1 $Kj 2, at each 
height (Ref. 2). 

In interpreting the numerical solutions of the present analysis, 
it should be pointed out that in addition to the assumptions previously 
discussed, equations (2) contain another source of potential error. 
Ideally, the unaccelerated flight of a single body can be completely 
simulated in the wind tunnel, whereas the separating flight of two 
parallel stages is impossible to simulate either kinematically or 
dynamically with fixed models in the test section. Thus, the direct use 
of the measured coefficients in equations (2) involves the error due to 
the wind-tunnel simulation technique where the static data are obtained 
for 7 ^ = 72 and (V c ) = (V c )^. This error is lessened by using the 
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instantaneous velocity of each stage in equations (3) ; however, even 
though the flight-path angle of each vehicle is varying during the 
separation maneuver the procedure does not take into account the effects 
of different flight-path angles on the static aerodynamic coefficients 
nor any scaling or Reynolds number effects . 


CHAPTER IV 


RESULTS AND DISCUSSION 

The results have been divided into two principal parts* These 
two parts consist of the experimental results obtained for the two stages 
in close proximity and the numerical results obtained by integrating the 
equations ot motion for each of the two stages during a staging maneuver. 
Since the aerodynamic interferences that occur when two bodies are in 
close proximity represent potential stability and control problems for 
both vehicle stages, the experimental results are used to discuss the 
implication of these interferences on the stability characteristics of 
the stages when in close proximity, while the numerical results are used 
to discuss the potential effects on the vehicle behavior of the aerody- 
namic int erf erence , 


Experimental Results 

The basic longitudinal aerodynamic characteristics of the first 
and second stages at different vertical and longitudinal spacings, 
incidence angles, and Mach numbers are presented in Figures 7 to 10. 
For reference, the interference-free data, or the data at a very large 


separation distance 




, are also shown on these Figures . 


First stage .- For the first stage at M = 3, i = 0 , and 

— = -0.051 (Fig. 7a) the static stability level, that is, the slope of 
l l 

the pitching-moment curve did not change appreciably with vertical 
spacing. However, there is a large positive increment in due to 
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the interference and a significant change in the angle for trim. Also, 

for the maximum vertical spacing tested, large interferences are still 

present because the curve at ^ = 0.224 has not approached the 

*1 

interference-free curve ^ = <c\. Figure 7a also shows that there are 

large negative increments in Cfj,l due to the interferences with no 
appreciable change in (C^) . Increasing i from 0° to 10° at M = 5 

and — = -0.051 (Figs. 7a, 7b, and 7c) caused a large positive incre- 

!l / 
ment in C m ^ and a large negative increment in C N,1 * Compare , 

' Az \ V 
N ^ at — = 0.127, for example. 1 

At the other longitudinal spacing, — = 0.l6o, and i = 0° 


and C 


(Fig. 7d) not only did the magnitude of C m ^ vary with vertical 
spacing, but also the stability level decreased with increasing vertical 
spacing such that at — = 0.199 and 0.224, the first stage is statically 

h 

unstable for less than about 2°. Figure 7d also shows that in 

comparison to the interference free curve, both (Cj^) and the 

magnitude of ^ a l a i greater than about -8 decreased due to the 

presence of the second stage. However, increasing the vertical spacing 

did not appreciably change either (Cjj ) or Cjj Increasing i at 
a a 1 

— = 0.160 (Figs. 7d, 7e, and 7f) caused further decreases in the 
Z 1 

stability level, magnitude of an d magnitude of Cj^i as compared 

to the interference-free curves. 

At M = 6, i = 0°, and ^ = -0.051 Figure 8a shows that both 
the magnitude of and stability level of the first stage varied 

with vertical spacing. This is in contrast to the data at M = 3 
(Fig. 7a) where only the magnitude of i was shown to vary with 
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vertical spacing. Also in contrast to the data at M = 5 the data at 

the largest vertical spacing tested, ^ = 0.224, at M = 6 (Fig. 8a) 

l l 

indicates that the first stage is approaching interference-free condi- 
tions since both the normal-force and pitching -moment curves are 
approaching the interference-free curves . However, at this vertical 

A *7 

spacing and also at — = 0.199> the first stage has pitch-up tendencies 

h 

at a-2 greater than 3°. Increasing i from 0° to 10° (Figs. 8a, 8b, 

and 8c) at M = 6 and — = -0.051 caused positive increments in 2 

l l * 

and a degradation of the static stability at all vertical spacings. 

Figures 8a, 8b, and 8c also show that the presence of the second stage 

caused negative increments in Cjj 2 with some localized changes in 

(Cw ) a.:, compared to interference-free conditions. 

1 AX 

At the other longitudinal spacing, — = 0.160, and M = 6 

h 

(Figs. 8d, 8e, and 8fi, the first stage was statically unstable at the 
smaller vertical spacings for <3.2 less than about 2°. These same 
figures also show that the magnitude of C N and (Cm) 

TX 2 

varied with vertical spacing. 

Second stage .- For the second stage at M = 3> i = 0°, and 
= -0.051 (Fig. 9 a), both the stability level and magnitude of (^2 

il > 

A 7 

varied with vertical spacing. At — less than 0.175 , the second stage 

l l 

was approximately neutrally stable or statically unstable at all 02 . 

Coupled with these changes in stability and magnitude of 0^2 are 

changes in both C„ Q and (C, T ) . Increasing i from 0° to 10° at 

5 

Ay ^ 

M = 3 and -r— = -0.051 (Figs. 9&, 9b, and 9c) caused the second-stage 
l l 

stability level to decrease and produced large negative increments in 
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the magnitude of C m ^ 2 as well as large positive increments in 
Comparison of these data (Figs. 9b, 9c) with the corresponding data 
for the first stage (Figs. 7a, 7b, and 7c) indicates that increasing i 
would tend to cause the first stage to nose up but away from the second 
stage, while the second stage would tend to nose down but away from the 
first stage. 

A ^ 

Increasing — from -0.051 to 0.160 for the second stage at 

h 

M = 5 (Figs. 9d, 9e, and 9f) caused the vehicle to become statically 
stable at almost all vertical spacings, incidence angles, and a^. 
However, the presence of the first stage caused increments in both 
C m 2 and % 2 with the increments becoming larger with increasing i. 

Comparison of the data for the second stage at M = 5 (Fig. 9) 
with the data at M = 6 (Fig. 10) at similar geometric conditions, 
indicates that the results at both Mach numbers are very similar. 

However, increasing the Mach number from 5 to 6 decreased the actual 
magnitude of the force and moment coefficients . 

Because of the large interference increments found for the second 
stage, an investigation was made to determine the control effectiveness 
of the second stage at interference-free conditions. Figure 11 presents 
this control effectiveness data for control deflections of 0° and ±25°. 
Comparison of this data with the data in Figures 9 and 10 indicates that 
there appears to be no reasonable control deflection which could overcome 
all of the large pitching -moment increments shown for the second stage 
when in proximity to the first stage. Furthermore, the control 
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effectiveness of the second stage may degrade due to the interferences 
or when in proximity to the first stage. 

Comments on interference increments .- Further analysis of the 
data shown in Figures 7 to 10, along with schlieren data, has led to 
two major conclusions as to the source o^ the behavior of the force and 
moment coefficients as compared to their interference-free values . These 
can be illustrated in Figure 12 . The first conclusion is that the 
effects on the first stage are caused by the impingement on the first 
stage of the bow shock wave generated by the second stage. Consequently, 
the changes in forces and moments previously shown for the first stage, 
are approximately proportional to the strength of this shock wave and 
the area affected. The second conclusion is that the effects on the 
second stage are caused by the flow field from the first stage to which 
must be added the effects of the first reflection of this bow shock 
wave. Furthermore, not only are the static aerodynamic data for both 
stages dependent on the relative position and attitude of the second 
stage and Mach number; but the data are also dependent on the relative 
sizes of the two stages, since the flow field of the first stage and the 
strength of the bow shock wave generated by the second stage are direct 
functions of the size of the two vehicles. The relative size of the 
two vehicles depends on the take-off mode of the launch vehicle (vertical 
versus horizontal), mission requirements, and the choice of propulsion 
systems for the launch vehicle. 
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Numerical Results 

When the equations of motion (eq. (l)) for each stage are coupled 
through the aerodynamic coefficients, the result is a system of 12 simul- 
taneous ordinary linear differential equations with 12 initial condi- 
tions. Since the initial time is conveniently chosen to be the moment 
of release, not all 12 initial conditions are independent. For instance, 
at release both stages are still flying together; consequently the 
initial velocities, both in magnitude and direction, must be the same. 

The initial spacing distance places similar requirements on the initial 
location and angular orientation of the two stages . 

The equations of motion were integrated numerically on a digital 

computer. The static aerodynamic inputs into the program consisted of 

the normal -force and pitching-moment coefficients for the two stages at 

different relative positions and attitudes when the stages were in close 

proximity to each other. The axial -force data used was that which was 

obtained at interference-free conditions for each stage. The dynamic 

derivatives C N , C , and CL. were estimated for each stage from 
fl q “q ®w 

data for similar type configurations such as the XB-70 for the first 
stage and the HL-10 for the second stage. However, as will be illustrated 
later, safe separation trajectories were not obtainable with these values 
and larger values had to be selected to achieve safe separation. Vehicle 
input data consisted of the masses of the individual stages, and 

m 2 , and the pitch moments of inertia, (iyy)^ and (lyy) 2 * Other 
vehicle input data consisted of the initial geometric quantities 
describing the initial position and attitude of the second stage with 
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respect to the first stage, namely, i. Ax, and Az. Quantities 
characteristic of the trajectory flown by the launch vehicle at staging 
conditions were also treat 3 d as inputs and included the altitude, h, 
pitch velocity, q, flight-path angle, 7 , velocity, V c , and angle of 
attack, a-^ . 

For the numerical integration the size of the computing interval 
was carefully selected. This selection was accomplished by the use of 
two different integration techniques for each of which the size of the 
computing interval was systematically varied. Comparison of the solu- 
tions obtained ir. this manner indicated that a At of 0.005 seconds 
resulted in a negligible error in the numerical results. 

The computer time involved in obtaining solutions puts practical 
limits on the number of variables and vehicle characteristics which can 
be investigated. Accordingly, for the present analysis, staging was 
assumed to occur in climbing flight; and Mach numbers of 3 and 6 were 
selected for the analysis. The Mach number 3 condition could correspond 
either to an abort or to a normal mission staging condition, while the 
Mach 6 condition would probably be a mission staging condition only. 
These considerations led to the selection of a number of fixed vehicle 
characteristics and initial conditions given in Table II. 

A typical altitude Mach number ascent trajectory plot, shown only 
for Mach numbers less than 8 , is illustrated in Figure 13 and indicates 
where abort and mission staging conditions might be required. Typical 
reusable launch vehicles employing air-breathing engines only would 
probably stage at the lower altitudes and resulting higher dynamic 
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pressures , whereas rocket -powered vehicles would stage at the higher 
altitudes and lower dynamic pressures. Normal mission staging could 
occur, depending on the vehicle concept, anywhere above a Mach number 
of about 3; while abort staging could occur at any Mach number below 
mission staging. Abort separation for either rocket or air-breathing 
type vehicles would likely occur at the higher dynamic pressures. 

Example solutions .- Figure l4 illustrates typical numerical 
results at M = 3 for two particular choices of the initial conditions 
and the damping in pitch characteristics. The vehicles were assumed to 
have safely separated when the vehicles were one first-stage body length 
apart and to have collided when the extremities of the vehicles touched 
each other. The main points of this Figure Eire that (l) the time 
required for the vehicles to collide or separate was of the ore e' of 
3 seconds, a very short time for aerodynamic controls to be effective 
(probably of the order of 2 seconds at these speeds and altitudes); 

(2) the net change in the initial and final velocities of either stage 
during the separation maneuver would amount to a net Mach number change 
of about 0.10 and would have a negligible effect on the static aerody- 
namic coefficients, thus Justifying assumption 1 in Chapter III; 

(3) the difference between y and 7 2 which was a maximum of 1.75° 
for the safe separation case, could introduce some errors in the calcu- 
lations (since the static aerodynamic coefficients were obtained for 

7 = 7^) and could affect whether the vehicles would collide or safely 


separate. 
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As has been previously implied, no attempt was made in the present 
study to investigate the effect of all important vehicle characteristics, 
initial conditions, and damping characteristics either separately or in 
all reasonable combinations. Instead, attention was given to some of the 
least accurately known damping variables, namely Qr id ( c m q ) 2 > 

and to some of the initial conditions which have obvious bearing on 
successful separation. Included in the latter category were the Initial 
relative attitude, dynamic pressure, and vertical thrusting. The deunping 
variables (C^) and (Cm^) are to obtain either analytically 

or experimentally at these high Mach numbers even for the interference- 
free condition. Since the static aerodynamic data indicated large changes 
in the pitching -moment characteristics for both vehicle stages when they 
were in close proximity, it can only oe expected that would also 

vary during the separation maneuver. The manner in which varies 

when the vehicles are in close proximity is beyond the scope of the 
present investigation and consequently, as was mentioned previously, 

C m ^ was held constant during the separation maneuver. 

Effect of dynamic derivatives .- Figure 15a shows the effect of the 
dynamic derivatives at M = 5 and h - 70,000 feet where (C^)^ is 
plotted against Indicated in zhis Figure are the regions where 

the two vehicles collided, an approximate region where they are safely 
separated, and a region of uncertainty as to whether the two vehicles 
would collide or separate. Also indicated in this Figure are the best 
current estimates of the interference-free values of for both the 

first and second stages . These values are one or more orders of 
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magnitude smaller than the values needed to achieve safe separation 
at these conditions . Figure 15b indicates similar results at M = 6 
and h = 110,000 feet. Although the magnitude of and for 

both the first and second stage have not been perturbed in this analysis, 
it can be expected that similar figures would be generated for these 
quantities at both Mach numbers. 

Effect of initial attitude .- For the purpose of illustrating so^e 
of the other important variables that may need to be considered to achieve 
safe separation, Figure lo shows the type of results that can be 
generated by selecting a value of which lies in the safe separation 

region for both the first and second stage at M = 5 and 6 (square 
symbol. Fig. 15). Here the incidence angle and the angle of attack of 
the first stage have been varied at M = 3 (Fig. l6a) and M = 6 
(Fig. 16b). Again, the regions of collision, uncertainty, and safe 
separation are shown. Similar results are again illustrated at both 
Mach numbers . 

Effect of dynamic pressure .- Figure 17 shows the effect of varying 
the dynamic pressure or altitude for the separation maneuver by taking 
a suitable safe separation value of i and a ± (circular symbol, 

Fig. 16) . The distance between the centers of gravity of the two stages 
in 2 seconds after release has been plotted against dynamic pressure. 

The dynamic pressure range is for altitudes above and below the nominal 
altitude for a rocket -powered vehicle at M = 5 and 6. Two curves are 
presented; one is for the condition in which no thrust was used in the 
calculations and the other is for the conditions in which a hypothetical 
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vertical downward thrust corresponding to a thrust -to-weight ratio of 
0.1 has been applied to the first stage. Using either vertical or 
longitudinal thrust would present plume-body interactions which, in these 
calculations, have not been evaluated. For the no-thrust condition at 
M = 3 (Fig. 17a) and zero dynamic pressure, the position and attitude 
of the vehicles do not change; and separation does not occur. Safe 
separation can be expected only at the higher dynamic pressures for the 
no-thrust condition. For the vertical -thrust condition at M = 3, safe 
separation could result at low dynamic pressures but as the dynamic 
pressure increases, there is a region of collision. At the higher 
dynamic pressures, safe separation would also be predicted. 

At M = 6 (Fig. 17b) safe separation was achieved for all 
conditions except at the extremely low dynamic pressures using no 
vertical thrust . However, the values of required for both the 

first and second stage to achieve safe separation are much larger than 
the best current estimates of the interference-free values. 

Figures 15, 16, and 17 have indicated some of the potential 
problems associated with parallel stage separation. These results have 
indicated that evaluating the potential of any parallel staged system 
will require a critical trajectory analysis of the two vehicles when in 
close proximity to determine the constraints under which safe separation 
of the two vehicles may be achievable. 


CHAPTER V 


CONCLUDING REMARKS 

An exploratory study has been made of the stage separation of 
parallel staged reusable launch vehicles . Static longitudinal aerody- 
namic data were obtained for both stages of a representative two-stage 
rocket -powered reusable launch vehicle concept when the stages were in 
close proximity to each other. The effects of vertical spacing, longi- 
tudinal spacing, and incidence angle were determined at Mach numbers of 
3 and 6. Using the wind-tunnel data together with estimates of the 
dynamic derivatives obtained for similar type vehicles and also the 
vehicle and initial trajectory input data, the implication of the inter- 
ferences on the individual vehicle dynamics was obtained at Mach numbers 
of 3 and 6 by numerically integrating the equations of motion for both 
vehicle stages . 

The experimental results indicated that for the first stage, the 
proximity of the second stage produced large positive increments in the 
pitching -moment coefficient coupled with large negative increments in 
the normal -force coefficient. Increasing the incidence angle of the 
second stage resulted in further positive increments in the pitching - 
moment coefficient and further negative increments in the normal -force 
coefficient. Ir general, the presence of the second stage caused the 
first-stage static stability level to decrease and at some of the 
geometric conditions tested, the first stage was statically unstable. 
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For the second stage, the proximity of the first stage generally 
resulted in large negative increments in the pitching coefficient and 
large positive increments in che normal -force coefficient. Increasing 
the incidence angle resulted in further negative increments in the 
pitching -moment coefficients and further positive increments in the 
normal -force coefficients. The presence of the first stage generally 
caused the second stage to become unstable at almost all geometric 
conditions tested. 

From control effectiveness data obtained for the second stage at 
interference-free conditions, there appears to be no aerodynamic control 
surface of reasonable dimensions which coula overcome the large pitching- 
moment increments . 

The observed changes in forces and moments were primarily caused 
for the first stage by the bow shock wave generated by the second stage. 
For the second stage, the changes in forces and moments were primarily 
caused by the flow field from the first stage to which must be added the 
effects of the first reflection of the bow shock wave generated by the 
second stage. 

The static aerodynamic data for both stages were found to be 
dependent on the relative position and attixude of the second stage, the 
flight Mach number, and the relative sizes of the two stages. The 
relative sizes of the two stages depend on the take-off mode of the 
launch vehicle (vertical versus horizontal), mission requirements, and 
the choice of propulsion system for the launch vehicle. 
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The numerical results indicated that safe separation is strongly 
dependent on the dynamic derivatives, the initial attitude of the 
two stages, and the dynamic pressure at staging. Perhaps the most 
significant result observed is that to achieve safe separation in the 
sensible atmosphere, the damping in pitch values required are one or 
more orders of magnitude larger than the best estimates for their 
interference-free values. However, to realistically assess the impact 
of the dynamic derivatives, dynamic derivative data are strongly needed 
for these types of vehicles at high supersonic and hypersonic Mach 
numbers; at interference-free conditions as well as when the vehicles 
are in close proximity to each other. 

These results have indicated that a potentially hazardous situa- 
tion could be expected in separating parallel -arranged stages. However, 
to evaluate the potential of safely separating any given system will 
require critical trajectory analysis of the system under consideration 
to determine the constraints for which safe separation may be achievable. 
Furthermore, other avenues of approach such as trajectory shaping, 
vehicle shaping, various thrusting maneuvers (including attitude control 
thrust), stability augmentation, and other auxiliary devices need to be 
considered. 

The experimental technique and method of analysis used in the 
present study are applicable to the separation of external stores from 
aircraft flying at high supersonic or hypersonic speeds. There are, 
however, two limitations which might present some difficulties. The 
first of these is that the experimental aerodynamic data cannot be 
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obtained, at present, if any lateral motion occurs during the separation 
maneuver; therefore, it was assumed that the separation of both vehicles 
takes place in the common plane of symmetry. The second of these limita- 
tions is that, even though the flight-path angle of each vehicle is 
varying during the separation maneuver, the procedure does not take 
into account the effects of different flight -path angles on the aero- 
dynamic coefficients; the aerodynamic coefficient data are obtained for 
the condition of flight -path angles of both vehicles equal to each other. 
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tab:j: t . - geometric characteristics of models 


AND REFERENCE DIMENSIONS 

First stage: 

Body length 23.760 inches 

Overall length ... 23.760 inches 

Aspect ratio 1 .865 

Span 13.802 inches 

Total wing area 102 .135 inches^ 

Reference area 102.135 inches 2 

Root chord . 14.800 inches 

Tip chord 0 

Mean aerodynamic chord 9 • 86 J inches 

Reference length 9-867 inches 

Second stage: 

Body length 11.465 inches 

Overall length 13.706 inches 

Aspect ratio O .678 

Span 6.050 inches 

Planform area less vertical fins 53-980 inches 2 

Reference area 53-980 inches 2 

Reference length 11.465 inches 
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TABLE II VEHICLE CHARACTERISTICS AND INITIAL CONDITIONS 


First stage: 



S . 

c n*w 

Craq 

% 

>1 


. . . 3^5,000 pounds 

11.00 x ICP slugs-ft 2 

78.11 feet 

6,401 feet 2 

0 

variable 

0 

26° 


Second stage: 




& 

*1 

Az 

Z 1 


?2 


. . 303 >000 pounds 

1.75 X 106 slugs-ft 2 
.... 90.76 feet 

.... 3,383 feet 2 

0 

variable 

0 

0.051 

0.127 

26° 
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APPENDIX 


DERIVATION OF THE RIGID-BODY EQUATIONS OF MOTION 
FOR AN AIRCRAFT WITH THE CONTROLS FIXED 

In the interest of completeness, the rigid-body equations are 
derived from first principles. The method is similar to the approach 
used by Etkin (Ref. 3)- 

Let 6m (Fig. A-l) be an element of mass of the airplane and 
have a velocity v relative to space-fixed axes. Let 5F be the 
resultant force vhich acts upon it. 


5F = 6m ^ 
at 


Summing all elements of the airplane gives 

f? “ at £ 7 ^ 


(A-l) 


(A-2) 


where 


_ _ ds 

V " V c + dt 


(A-3) 


First, look at the term ) v 6m 


I 5 


) v 60 = ) (v c + |j)6m - "So + ^ I 


s 6m 
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Since the equation? are being written for a coordinate system at the 
mass center, \ s 6n: - 0. Therefore, 


^ v 5m = mv,, (A-4) 


Substituting equation (A-*t) into (A-2) gives 


I 



or 


T = 



(A-5) 


It is now desired to obtain an equation for the moments about the mass 
cen + er. The moment of momentum of 6m is, by definition, 

6H = s x v 6m 


Consider 


-^-(6H) = -^-(s x v)5m = ^xv6m + 6x^6m 
dt du dt dt 


(A-6) 


From equation (A-5) 


ds 

dt 


v 


- v c 
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Also. 


sx?r8m=8x5F=6G 

dt 


Therefore, equation (A-6) now becomes 


6G = t"( 5H) - (v - v c ) x v 6m 
dt 


or 


6G = t"( 8H) + v c x v 6m 
dt 


(A-7) 


Equation (A-7) is now summed for all elements 


)> • S > 


(6H) + X ) V6n = ^ 


(51) + V x mV c 


or 


»-£ > 


(A-8) 


Now, \ 6G = the resultant external moment about center of gravity 

while ) 6H -- the angular momentum of the airplane about center of 
gravity. Consequently, equation (A-8) becomes 


p _ dH 
G " dt 


(A-9) 



Ix is necessary now to evaluate the angular momentum H. 


H = ^ 5H = ) (sx v)5m 

Let the angular velocity of the airplane be 

g 5 = pi + qj 4 rk 

The derivative of a vector E, referred to a frame of reference 
rotating with angular velocity o> can be shown to be given by the 
following equation: 


dE 

dt 



Moving system 


+ CJD x E 


Consequently, the velocity of a point in a rotating rigid body is 
given by 


as 

dt 


df 

dt 


(Moving system 


+ 03 x s 


or 


V = V c + (D x s 


Therefore, 


(A-10) 


(A-ll) 


(A-12) 
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H = ) s x v„ 5m + ) s x (cjo x s)5m 


Now, since ) s 6m = 0, the expression 


^ s x v c 5m = ^ s 6m x v c = 0 


The second term on the right-hand side of equation (A-13) becomes 


^ s x (<cx s)6m = ^ jjo(s 2 ) - s(o> • sf]6m 


Consequently, equation (A-13) reduces to 


H = a) s 2 6m - ^ s(o) ♦ s)6m 


With s = xi + yj + zk and H = H x i + Hyj + H z k, equation (A-l4) 
broken down into its scalar components which are as follows : 


H, 


= p ^ (x 2 + y 2 + z 2 )6m - ^ (px 2 + qxy + rxz)5m 


Hy = q ) (x 2 + y 2 + z 2 )6m - (pxy + qy 2 + rys 


)6m 


H z = r 


^ (x 2 + y 2 + z 2 )6m - ^ (pxy + qyz + rz 2 )5m 


(A-13) 


(A-14) 
can be 
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or 


h x = p N (y I 2 + z 2 ) om - \ (qxy + rxz)6m 


Hy = q J (x 2 + z 2 )6m - ) (pxy + ryz)6m 


H z = r \ (x 2 + y^m - ) (pxz + qyz)6m 


Nov 


i xx. = ) (y 2 + z2 ) 6m 


Iyy s y (x 2 + z 2 )5m 
I zz = ) i* 2 + y 2 ) 5m 


I 


xy 


xy 6m 


I = \ xz 6m 

xz /_, 


I 


yz - 


)_ 


(A-15) 


yz 6m 
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Therefore, equations (A-15) can be written as 

H x = x xxP - ^ ' X xz r 

H y = “ X xyP + X yy q " X yz r (A-l6) 

H z = - X xzP - Iy.q + X zz r 

The vector equations of motion (eqs. (A-5) and (A-9)) then become, when 
referred to the coordinate system fixed to the airplane (X, Y,Z) : 



Moving system 


+ 


m(o> x v c ) 


(A-17) 


+ <5 X H (A-l8) 

Moving system 

where 

F = F x i + F y j + F z k 
G = M X I + M y ' j + M z k 
H = H X I + H y j + H z k 
o> = pi + q J + rk 
v c = ui + vj + wk 

Equations (A-17) and (A-l8) have the scalar components: 
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F x = m(u + qw - rv) 

Fy = m(v + ru - pw) (A-19) 

F z = m(w + pv - qu) 

M x = S x + qH z - rHy 

My = Hy + rH x - pH z (A-20) 

M z = H z + P Hy - qH x 

Since the frame of reference for the equations of motion is fixed 
to the airplane and moves with it, the position and orientation of the 
airplane cannot be described relative to it. For this purpose Euler's 
Angles are used to describe motion of the airplane. The Euler ian angles 
are defined as the three successive angles of rotation needed to carry 
out the transformation from a given cartesian coordinate system to 
another. The following three rotations for the Euler's Angles will be 
used (see Fig. A-2): 

(a) A rotation i|r about Z^-axis. 

(b) A rotation 9 about Y 2 -&xis. 

(c) A rotation 0 about X^-axis. 

To obtain the necessary transformation matrix, a rotation about Z^- 

axis is done first. The transformation matrix is as follows: 





Figure A-2 


A = 


cos \|r -sin i|r 0 

sin 'Jr cos \|r 0 

0 0 1 


Secondly, rotate about the Yg-axis . The transformation matrix is as 


follows : 
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B 


cos 0 0 sin 0 

0 10 

^-sin 0 0 cos 0^ 


Thirdly, rotate about the X^-axis. The transformation matrix is as 
follows : 


C 


10 0 
0 cos 0 -sin 0 

_0 sin 0 cos 0 


The product matrix D = ABC then is as follows, where D = the trans- 
formation matrix from body coordinates to space coordinates. 


D = ABC = 


cos \| ( cos 0 sin 0 sin 0 cos a|< 

- sin cos 0 

+sin \|( cos 0 +sin sin 0 sin 0 
+ cos \|> cos 0 


-sin 0 


sin 0 cos 0 




* — 

+ sin ^ sin 0 

+sin \|r cos 0 sin 0 
- sin 0 cos \|< 

cos 0 cos 0 


(A-21) 


Consequently, to obtain the transformation from space coordinates to 
body coordinates, we use the transformation matrix D 



where jjcQ refers to the space coordinate system and QTJ refers to the 
body coordinate system. To obtain the coordinates of the flight path 


4i 


relative to a fixed frame of reference as shown in Figure A-2, use the 
transformation matrix D. 

If the velocity matrix in the fixed frame is 


u’ 


dx' 

dt 

v-i 


dy' 

V 


dt 

v' 


dz' 

dt 


the relationship between the body coordinate system and fixed coordinate 
system is found as follows : 


_ dx' 


u 

dt 


5Zl 

dt 

■H 

V 

dz' 


w 

dt _ 




where u, v, w are the components of the velocity vector in the body 


coordinate system, or 



k2 


= u cos 9 cos \(i + v(sin 9 3in 0 cos \|i - sin ^ cos 0) 
dt 

+ w(cos \|i cos 0 sin 9 + sin >|r sin 0) 

= u sin cos 9 + v(sin \|/ sin 0 sin 9 + cos cos 0) v 
dt / 

+ w(sin f cos 0 sin 9 - sin 0 cos i|r) 

= -u sin 9 + v sin 0 cos 9 + v cos 9 cos 0 

dt j 


(A-22) 


It is now desired to express the orientation of the airplane in 

• • r 

terns of the angular velocity components (i|r,9,0). From the rotations 
used previously, we have 

= 0 = angular velocity about X-axis 

= 9 = angular velocity about Y 2 -axis (A-2J) 

o> * = angular velocity about Z^-axis 

In vector notation, 

• • • 

o> = 10 + J^e + k p v 

In the body axes coordinate system (X,Y,z) (see Fig. A-2), 

o> = pi + qj + rk 

If the components of given by equation (A-25) are projected onto the 
X,Y,Z axes, or body axes, the following relations are obtained in matrix 


notat ion 
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^^body 


((JO ) 

V 0 'body 


1 

■“O" 


_ -ijf sin 0 


fr]' 1 

0 

= 

sin 0 cos 

0 




Jjr cos 0 cos 

0. 



“0” 


“0 

- 000 

0 

= 

0 cos 0 


^0 _ 


^-0 sin 0_ 


$ 


(“#) 


body 



Consequently, 



■ K’bod, + ‘^'body * ( Vbody 


or 


p = 0 - sin 0 + 
q = 0 cos 0 + ^ sin 0 cos 0 
r = cos 0 cos 0-0 sin 0 


(A-24) 


For the problem at hand, the body axis system was chosen, in 
which case the X-axis is fixed to a longitudinal reference line in 
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the airplane. The airplane was assumed to have exact symmetry about the 
selected axis system. In this case, 

x vz “ Ixy = 0 (A-25) 

With the aerodynamic forces (including propulsive forces) denoted 

by F', F', F', the resultant external forces are (see Fig. A-5): 

x y z 

F = F ' - mg sin 9 

X A 

Fy = Fy + mg cos 0 sin 0 (A-26) 

F z = F^ + mg cos 6 cos 0 

The kinematical and dynamical equations derived in the foregoing 
are now collected. 

Substituting equations (A-19) into (A-26) and solving for u, v, 
w, the following equations are obtained 

u = rv - qw - g sin 0 + j Fj 

v = pw - ru + g cos 9 sin 0 + - F' (A-27) 

my 

w = qu - pv + g cos 9 cos 0 + ; f 

m z 

Substituting equations (A-l6) into (A-20) and solving for p.. q, 
r, the following equations are obtained 
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where 


P = (A^pq - (Ag)qr + (Aj)^ + (A^)^ 
q = (A^)pr - (AgHr 2 - p^) + (Ay)M y (A-28) 

r = (Ag)pq - (AjJqr + (A^jR^ + (A^)^ 


b ■ 


I I 

- I 2 



XX zz 

xz 



T (J 

- I 

+ 

I ) 

xz ' xx yy 


zz y 

I 2 - 

I I 

+ 

H 

no 

XX 

xx yy 


xy 

1 

CVJ 

h- 

I I 

+ 

I 2 

zz 

zz yy 


XZ 

u zz - 

l-xx) 




L yy 


a 6 

\ 

a 8 

*9 



The equations used to describe the angular positions of the airplane 
are equations (A-24) - solved for 6, 0, and 4 -: 



^7 


9 = q cos 0 - r sin 0 

0 = p + q sin 0 tan 9 + r cos 0 tan 9 (A-2?) 

\> = (q sin 0 + r cos 0)sec 9 

The position of the center of gravity of the airplane referenced 
in earth-fixed coordinate system is given by equations (A-22): 

"V 

x' = u cos 9 cos + v(sin 0 sin 9 cos - sin i|r cos 0) 

+ w(cos t cos 0 sin 9 + sin \)f sin 0) 

y' = u sin ^ cos 9 + v(sin sin 0 sin 9 + cos ^ cos 0) \ (A-30) 

+ v(sin + cos 0 sin 9 - sin 0 cos 

z' = u sin 9 + v sin 0 cos 9 + w cos 9 cos 0 



^ o 
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Figure 1.- Schematic of 2-foot hypersonic facility testing section. 



Figure 2.- Details of models. All dimensions are in inches 





Figure 2 Continued. 
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Figure 2.- Concluded 








54 



Second -stage 
release 
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Figure 6.- Staging sequence. 






Figure 7 »- Continued . 
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(d) i - < 3 °j ~ • 0.160. 


Figure 8.- Continued. 






(e) i ■ 5°| “ - 0.160. 
l l 


Figure 8.- Continued. 
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Figure 9-- Longitudinal aerodynamic characteristics of the second stage 

at a Mach number of 3. 
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(d) 1-0°; ^ - 0.160. 

*1 

Figure 9*- Continued. 
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Figure 9*- Continued 
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Figure 9** Concluded. 
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(a) i - 0°# “ -0.051. 

*1 

Figure 10.- Longitudinal aerodynamic characteristics of the second stage 

at a Mach number of 6. 
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(d) 1 ■ 0°; y- ■ 0 . 160 . 

Figure 10.- Continued. 
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(a) M = 3. 

Figure 11.- Control effectiveness of the second stage at interference 

free conditions . 
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(b) M = 6. 


Figure 11.- Concluded . 




Figure 13.- Staging conditions. 
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# Estimated interference -free values 
■ Point for further analysis 



(a) M - 3; h - 70,000 ft; a x » 6°; 1 - 5°. 

Figure 15.- Effect of dynamic derivativea on the separation maneuver. 
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~V. 


• Estimated interference -free values 
■ Point for further analysis 



0 -20 -40 -60 

( c m q ), ,rad~* 


(b) M ■ 6j h - 110,000 ftj ■ 6°| i * 5°« 
Figure 1?.- Concluded. 
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O Point for furthCi analysis 



(b) M * 6# h - 110,000 ft# (Cfoq)^ » -^0 rad" 1 # (0^)^ = -5 rad" 1 . 

Figure 16 .- Concluded. 


Vertical thrust 



(a) M = 3j ai = 6°; i = 5°; (C^ = -4o rad’ 1 ,- (C^ = -20 rad 

Figure IT.- Effect of dynamic pressure on the separation maneuver. 
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Vertical thrust 



Figure 17.- Concluded. 


(b) M - 6f 
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